Evidence concerning Drying behavior of Ne near the Cs surface 
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Using density functional (DF) and Monte Carlo methods, we have studied the properties of Ne 
adsorbed on a Cs surface, focusing on the region at and near saturated vapor pressure (SVP). In the 
case of Ne/Rb, the experimental data of Hess, Sabatini and Chan are consistent with the calculations 
based on an ab initio fluid-substrate potential, while in the Ne/Cs case there is indication that the 
potential is ~ 9% too deep. In that case, the calculations yield partial drying behavior consistent 
with the experimental finding of depressed fluid density near the surface, above SVP. However, we 
find no evidence of a drying transition, a result consistent with a mean fleld calculation of Ebner 
and Saam. 

PACS numbers: 68.10.-m , 68.45.-v , 68.45. Gd , 68.35. Rh 



The subject of wetting has stimulated both basic and 
applied research for nearly two centuries. The Young 
equation implies the possibility of several alternative va- 
rieties of wetting behavior as a function of the adsorption 
system and the temperature (T). For example, at satu- 
rated vapor pressure (SVP), a liquid drop may bead up 
on a surface exhibiting a finite contact angle Q. This 
" non- wetting" behavior occurs when the following rela- 
tion between three surface tensions is satisfied: 

Act| = \asv — o-si \ < o'lv (0 < 8 < TT, nonwetting) 

(1) 

Here v, I, and s refer to vapor, liquid and solid, respec- 
tively. An alternative possibility is " wetting" (sometimes 
called " complete wetting" ) 

Act > aiy (9 = 0, wetting) (2) 

in which case a very thick liquid film (in equilibrium) is 
spread uniformly across the surface. The transition be- 
tween these regimes occurs at the wetting temperature 
Tw- Such a transition, anticipated by general arguments 
by Cahn Q and Ebner and Saam [Q), has been seen in 
a number of instances where the gas-surface attraction 
is weak. The theory envisions a third possible scenario, 
" drying" , in the case of a very weakly attractive interac- 
tion. The criterion is: 

A(T < — (Tjt, (8 = TT, drying) (3) 

Such behavior is manifested, in principle, as the pres- 
ence above SVP of a thick region of vapor intervening 
between the surface and the asymptotic bulk liquid. A 
recent experiment of Hess, Sabatini and Chan ||^ (HSC) 
found some evidence for such drying behavior for Ne 
adsorbed on the Cs surface. This system is a promis- 
ing candidate for drying since the theoretical well depth 



D 25 K of the gas-surface interaction is signifi- 
cantly less than the well depth (e — 33.9 -fC) of the Ne- 
Ne pair potential. In this paper, we report calculations 
relevant to the drying behavior observed in the HSC ex- 
periment. We have employed two methods, a Density 
Functional (DF) approach and grand canonical Monte 
Carlo (GCMC) simulations; we have improved somewhat 
upon techniques used in our previous studies of wetting 
transitions [^^. The DF method involves an empirical 
free-energy functional written in terms of the density p{r) 
of the fluid as: 

F[p] = FhsIp] + P{r^p{r')ua{\r- r'\)drdr' 

+ j p{r)Vs{r)dr + lFu (4) 

Here Fhs is the free-energy functional for an inhomo- 
geneous hard-sphere reference system, the second term 
is the usual mean-field approximation for the attractive 
part of the fluid- fluid intermolecular potential Ua, while 
Vs{r) is the external, static adsorption potential due to 
the surface. Fid is the ideal (classical) gas contribution. 
For Fhs we use the non-local functional of Ref. 0. 

The functional (^ involves three parameters: (i) a cor- 
rection to the the HS diameter, (ii) a parameter enhanc- 
ing the effective potential well depth of the fluid-fluid 
interaction, and (iii) the coefficient 7 in the last term of 
(^. A previous version of the functional (^), involving 
only the two parameters (i) and (ii), has been applied 
with success to the study of the wetting properties of 
Ar and Ne on different surfaces In particular a re- 

markably good agreement has been found § between DF 
calculations for the Ne/Rb system and the experimental 
results both showing a first-order wetting transition 
between 43 and 44 K. The small revision used here is to 
let the coefficient 7 in Eq.(U) be a free parameter; the 
best fit values are very close, for each value of T, to the 
value 7 = 1 assumed previously. 
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These three adjustable parameters are fit to reproduce 
properties of bulk Ne. In particular, we require that the 
experimental pressure and densities for the homogeneous 
system are reproduced at liquid-vapor coexistence (see 
Fig. ^ upper panel), and moreover that the stability 
condition /i„ ~ involving the chemical potentials of 
the two phases, is also satisfied. 
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FIG. 1. Upper panel: Pressure-density phase diagram of 
bulk Ne at SVP, at various temperatures. The squares are 
experimental points at coexistence, which are reproduced by 
construction by our functional. Lower panel: Liquid- vapor Ne 
surface tension ct;„. The squares are experimental points, the 
triangles are our calculated values, and the dashed line is a fit 
to the experimental data which goes to zero as (1 — T/Tc)^ , 
with the exponent /3 = 1.27 given by Renormalization Group 
Theory. 



We then use the functional to study the Ne liquid- 
vapor free interface and confirm that it yields good agree- 
ment with the (known) Ne surface tension aiv (see Fig. 
|l|, lower panel). 

When studying fluid-surface interactions, it is impor- 
tant to employ an accurate adsorption potential. Ab ini- 
tio potentials for many adsorption systems were derived 
recently by Chizmeshya, Cole and Zaremba jJ](CCZ). 
Overall, these have proven to be reliable for predicting 
contact angle and wetting temperatures of He and Ne on 
several surfaces p|,p|,p^ . Using the Ne/Cs potential, our 
previous GCMC study [|j found negligible adsorption of 
Ne over the entire range (T < 43 K) of that study, that 
for which the correlation length ^ of critical fluctuations 
is smaller than the smallest dimension (29 A) of the peri- 
odically replicated simulation cell; this result is consistent 
with the HSC data for this system. To more closely ap- 
proach the critical point {Tc — 44.4 if), we have increased 
the cell dimensions to lOOAx 100Ax200A. Our previous 
DF study j|] found a wetting transition for Ne/Rb at 
Tw ~ 43 K, consistent with the HSC data. With the re- 
vised DF method described above we confirm the Ne/Rb 
result and find a wetting transition for Ne/Cs at similar 



temperature {Tw ~ 43 — 44 K) . This is seen in Fig. ^ 
as a crossing of the two curves AtT(T) (curve (a)) and 
criv{T) (for the details of the method used to compute 
the surface tensions aij — s,l,v), see Ref. |^). Such 
a prediction disagrees with the HSC findings for Ne/Cs, 
of nonwetting at all T. We conjecture that this discrep- 
ancy occurs because the CCZ potential is somewhat too 
attractive. We therefore modify the original Ne/Cs po- 
tential in such a way that the modified potential has a 
slightly smaller well depth than the original one, as can 
be seen in Fig. |^. When using the modified CCZ po- 
tential, the wetting transition disappears, as can be seen 
from the curve (b) in Fig. |^ which no longer crosses 
aiv{T), i.e. a rather small (~ 9%) correction to the the- 
oretical well depth brings the wetting behavior into con- 
sistency with experiment. Such an "error" in the CCZ 
potential is compatible with the uncertainties present in 
its derivation (such as the jellium model of the surface 
and an empirical damping procedure applied to the dis- 
persion part of the attraction). 
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FIG. 2. Calculated surface tensions for Ne on different sub- 
strates. The two solid lines show the liquid-vapor surface ten- 
sion, aiv Dashed lines labeled (a) and (b): Aa for the Ne/Cs 
system calculated with the CCZ potential (curve (a)) and the 
modified-CCZ potential (curve (b)), respectively. The other 
lines show Act for Ne on ultraweak 9-3 surfaces, with decreas- 
ing values of the potential well depth D, as described in the 
text. 

Values of the contact angle cos(8) = Aa/aiv can 
be immediately extracted from Fig. 0. For the mod- 
ified CCZ potential we find, for instance, O = 94° at 
T = 38if, 6 = 84° at T = iOK, and 6 = 64° at 
T — A2 K. Measurements of these angles should provide 
a direct test of our calculations. 

Fig. H presents results for Aa(T) even for quite small 
hypothetical well depths. To simulate such ultraweak 
(UW) adsorbing surfaces, we use for simplicity a tunable 
3-9 model potential 
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where C3 is kept fixed to the Ne/Cs value, while the well 
depth D is arbitrarily varied. Fig. || displays results for 
the dependence of Acr(r) on the well depth D for such 
UW surfaces. Note that none of these curves crosses 
the curve — cr;^. This means that no drying transition 
occurs below T,. even for the weakest interaction consid- 
ered {D ^ Q.hK), a conclusion that is consistent with a 
general (mean field) result of Ebner and Saam. That 
argument implies that a drying transition can occur only 
at Tc as a consequence of the long-range van der Waals 
attraction. In order to check our DF results, we per- 
formed extensive GCMC simulations and confirmed the 
absence of a drying transition. What is thus responsi- 
ble for the apparent drying behavior for Ne/Cs reported 
by HSC? Fig. ^ presents a set of density profiles, for 
the two systems investigated in HSC, Ne/Au and Ne/Cs, 
for different values of the chemical potential ^ both be- 
low and above the value fiQ at liquid-vapor coexistence. 
Values of fi < fiQ result in profiles that have the vapor 
density as the asymptotic density p{oo) (far from the sur- 
face) whereas when /i > /io the liquid density is reached 
far from the surface. We note that these DF results are 
quite consistent with the GCMC data simulated with the 
same conditions. From the profiles shown in Fig. ^ one 
can compute the surface excess coverage per unit area 



AN/A 



dz[p{z) - p(oo)] 



(here A is the surface area and the lower limit zq is taken 
at the zero energy turning point of the adsorption poten- 
tial). 
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FIG. 3. Ne/Cs interaction potentials. The solid line shows 
the original CCZ potential, the dashed line shows the modi- 
fied CCZ potential described in the text. 



Results for An = AN/A are shown in Fig. ||, both 
below and above SVP. Note that the excess coverage is 
positive below po even if the attraction is very weak, as 
in the case of Ne/Cs. Above /io, the excess coverage for 
Ne/Cs is negative, i.e. there is a deficiency in mass be- 
cause the attraction is so weak. A positive An is found 
instead for Ne/Au, due to the compression of the liquid 
layers close to the surface induced by the strong adsorp- 
tion interaction. The small excess below SVP is in accord 
with HSC data. The defect coverage for Ne/Cs above po 
is An ^ 2.2 x 10~^ g/cm^, corresponding to ~ 0.9ML. 
The expected shift in frequency in a quartz microbalance 
experiment is A/ = {Ap/nR)An ~ A Hz, i.e. at 
least one order of magnitude smaller than the deficit re- 
ported by HSC above SVP. We believe that this apparent 
"discrepancy" is misleading. The experiment measures 
missing mass by its effect on the resonant frequency shift 
of the shear wave of the microbalance. However, a con- 
tributing factor is that a fluid which interacts weakly with 
a surface exhibits significant slip, as reported by Thomp- 
son and Robins . Hence the reported mass deficiency 
implicitly includes the reduction in mass dragged by the 
oscillator (as well as the reduced density of the neighbor- 
ing fluid). 
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FIG. 4. Equilibrium calculated density profiles for Ne/Au 
and Ne/Cs at T = 43 K, respectively, for difi'erent values of 
the chemical potential both below and above coexistence, z 
is the coordinate normal to the surface, measured from the 
surface plane position. 



Wetting and drying phenomena in planar geometries 
have been extensively studied in the last decade by us- 
ing MD computer simulations and DF calculations. In 
particular, model system of LJ fluids interacting with a 
9-3 wall were commonly considered. A drying transition 
between the triple-point Tt and the critical temperature 
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Tc is invariably observed in these simulations for a suffi- 
ciently low value of the ratio between the fluid-substrate 
and fluid-fluid interaction strengths, although with some 
controversy concerning the order of this transition (see 
Rcf. Jl^ , p^ , and references cited therein). However, de- 
spite the long-range character of the fluid-wall poten- 
tial used in these simulations, the interaction is usually 
cut-off, for computational reasons, at a short distance 
(tc ^ {2 ~ 3) <7lj), thus resulting in an effective finite 
range fluid- wall interaction. Inclusion of the full range of 
such interaction would results, as clearly shown by our 
results, in the absence of such drying transition for any 

T<Tr. 
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FIG. 5. Surface excess coverage AN/ A close to coexistence, 
for Ne/Au and Ne/Cs. /lo is the chemical potential at SVP. 

In summary, DF and GCMC calculations flnd a wet- 
ting transition for Ne/Cs if the ab initio potential is used 
to describe the fluid-substrate interaction, in disagree- 
ment with experiments. A plausible, small change in the 
attractive potential well depth changes the prediction to 
nonwetting for all T. The resulting behavior of the mass 
excess per unit area An is qualitatively consistent with 
the experimental data. Below SVP, An is small and pos- 
itive. Above SVP, an apparent discrepancy in the mag- 
nitude of An is attributed to the fact that significant slip 
occurs due to the weak interaction and resulting low den- 
sity is suggested by molecular dynamics results of 



Finally, our results are consistent with a general predic- 
tion of Ebner and Saam that a drying transition is shifted 
to Tc by the long-range gas-surface attraction. Our re- 
sults (to be reported in detail elsewhere) should be tested 
by contact angle and ellipsometric measurements of film 
thickness for Nc/Cs. 
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